home
***
CD-ROM
|
disk
|
FTP
|
other
***
search
/
EnigmA Amiga Run 1999 March
/
EnigmA AMIGA RUN 35 (1999)(G.R. Edizioni)(IT)[!][issue 1999-03].iso
/
earcd
/
-archivi
/
-recent2
/
amhelios.lha
/
AmHelios
/
hc_delta.cpp
< prev
next >
Wrap
C/C++ Source or Header
|
1997-07-13
|
2KB
|
84 lines
////////////////////////////////////////////////////////////
//
// HC_DELTA.CPP - Hemicube Delta Form Factor Class
//
// Version: 1.03A
//
// History: 94/08/23 - Version 1.00A release.
// 94/12/16 - Version 1.01A release.
// 95/02/05 - Version 1.02A release.
// 95/07/21 - Version 1.02B release.
// 96/02/14 - Version 1.02C release.
// 96/04/01 - Version 1.03A release.
//
// Compilers: Microsoft Visual C/C++ Professional V1.5
// Borland C++ Version 4.5
//
// Author: Ian Ashdown, P.Eng.
// byHeart Software Limited
// 620 Ballantree Road
// West Vancouver, B.C.
// Canada V7S 1W3
// Tel. (604) 922-6148
// Fax. (604) 987-7621
//
// Copyright 1994-1996 byHeart Software Limited
//
// The following source code has been derived from:
//
// Ashdown, I. 1994. Radiosity: A Programmer's
// Perspective. New York, NY: John Wiley & Sons.
//
// It may be freely copied, redistributed, and/or modified
// for personal use ONLY, as long as the copyright notice
// is included with all source code files.
//
////////////////////////////////////////////////////////////
#include "hc_delta.h"
// Static delta form factor arrays
float HemiDelta::side_array[HC_ArrayDim][HC_ArrayDim];
float HemiDelta::top_array[HC_ArrayDim][HC_ArrayDim];
HemiDelta::HemiDelta() // Class constructor
{
int i, j; // Loop indices
double da; // Cell area
double dx, dy, dz; // Cell dimensions
double r, x, y, z; // Cell co-ordinates
// Initialize cell dimensions and area
dx = dy = dz = 2.0 / (float) FF_ArrayRes;
da = 4.0 / ((float) FF_ArrayRes * (float) FF_ArrayRes);
// Calculate top face delta form factors
x = dx / 2.0;
for (i = 0; i < HC_ArrayDim; i++)
{
y = dy / 2.0;
for (j = 0; j < HC_ArrayDim; j++)
{
r = x * x + y * y + 1.0;
top_array[i][j] = (float) (da / (PI * r * r));
y += dy;
}
x += dx;
}
// Calculate side face delta form factors
x = dx / 2.0;
for (i = 0; i < HC_ArrayDim; i++)
{
z = dz / 2.0;
for (j = 0; j < HC_ArrayDim; j++)
{
r = x * x + z * z + 1.0;
side_array[i][j] = (float) (z * da / (PI * r * r));
z += dy;
}
x += dx;
}
}